Vibrational spectroscopy of polar molecules with superradiance 
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We investigate cooperative phenomena and superradiance for vibrational transitions in polar 
molecule spectroscopy when a high optical-depth (OD) sample is studied. Such cooperativity comes 
from the build-up of inter-particle coherence through dipole-dipole interactions and leads to the 
speed-up of decay process. We compare our calculation to recent work [Deiglmayr et al., Eur. Phys. 
J. D 65, 99 (2011)] and find very good agreement, suggesting that superradiant effects need to be 
included in a wide variety of ultracold molecule setups including vibrational and rotational states. 



I. INTRODUCTION 

Ultracold molecules have recently become an intensely 
explored area due to their fascinating properties. Many 
proposals have been brought up, for example, to im- 
prove the precision measurement of fundamental con- 
stants P, or test the fundamental law of physics, e.g., 
Q. As expected, the question of how to gain better ac- 
cess to internal degrees of a molecule to help better ma- 
nipulation becomes an important step. The most com- 
monly used tool of analysis in molecule systems is spec- 
troscopy. However, for molecules the resolution of com- 
plicated spectral lines with necessary accuracy is crucial. 
Typically, the spectroscopic model for internal degrees 
of a molecule has also been applied to an ensemble, by 
assuming that the overall radiative behavior does not de- 
viate dramatically from the direct sum of individual re- 
sponses of the ingredient molecules. Even in the case 
of amplified spontaneous emission (ASE, also called "su- 
perluminescence") [1, Q, where the medium nonlinear- 
ity plays an important role, the inter-particle coherence 
in terms of genuine quantum manybody nature is usu- 
ally overlooked. Such quantum manybody coherence can 
modify the emission rate significantly, resulting in "su- 
perradiance", which was first predicted by Dicke in 1954 
Since then plenty of efforts have been made to inves- 
tigate such phenomena both theoretically 0, and 
experimentally [1, [l3 - [20| . The interest in superradiance 
has been revived especially during recent years due to 
the advances of controllin g a tomic, molecular, and other 
quantum optical systems |2l| - [2a |. 

Looking at the simplest model, Dicke considered a 
case where the ensemble of excited two-level particles 
is confined in space to a size much smaller than the 
transition wavelength. In the limiting case where we 
are allowed to regard the whole ensemble as a point- 
like system, the system can be described by a collec- 
tive spin J = N/2 with N the number of particles and 
the factor i is for two-level atoms analogous to spin- 
i particles. As a result the overall decay rate is now 
given by T{J,M) = ^ {J M\D+ \J M) , where 7 is the 
free-space spontaneous decay rate, the collective rais- 
ing operator Z?+ = |e)i(g| with particle index i and 
the lowering operator = (_D+)'I'; M is the index 



of the Dicke ladder with —J<M<J. This factor 
{JM\D+D~\JM} = (J - A/)(J + M - 1) quantifies the 
enhancement of the decay rate, contributing to the origin 
of superradiance d, [l^ . 

Note that in the above discussion for Dicke superradi- 
ance, the sample size is considered infinitesimally small 
compared to the transition wavelength A. In reality, the 
particle ensemble must have a finite spread in space. 
Also, the distance-dependent dipole-dipole interaction 
breaks the Dicke symmetry and causes considerable de- 
phasing. A quantity called cooperativity can be defined 
to characterize the competition of these two length scales: 
C = A/'A'^/(47r), where M is the volume density of parti- 
cles. As one can expect, superradiance is more obvious 
when C ^ 1 but getting suppressed when C decreases. 
However, according to further analysis 0, [134111 j is 
more accurate to characterize superradiance behavior by 
the optical depth OD = AfX^d with sample diameter d, 
consistent with the known fact that the superradiance 
can be directional and hence depends on the geometry of 
the sample. Recent experimental technology has allowed 
to realize high OD for atomic and molecular ensembles. 
Therefore, cooperative phenomena such as superradiance 
should be observable in such systems. Indeed, a study 
on rubidium Rydberg atoms dealing with Rydberg state 
relaxation has showed that the associated life time is con- 
siderably shortened given OD w 10^ [H, Is^l- Its decay 
behavior agrees quite well with the theoretical calcula- 
tion in our previous work [30] , suggesting the emergence 
of superradiance. 

In this paper, we want to discuss the same effects for 
polar molecules specially for vibrational spectroscopy for 
which the parameter regime is accessible by current tech- 
nology. We have developed an effective two-body for- 
malism that accounts for the actual dipole-dipole inter- 
action for finite OD, for which the Dicke picture fails to 
be valid. We also compare our calculation to a recent ex- 
periment on the vibrational relaxation of LiCs molecules 
jSlj . Even though in this experiment OD is rather small 
such that the superradiance feature is not dominant, we 
observe a discrepancy between the measured data and 
the results given by the traditional single-particle rate 
equation approach including blackbody radiation. Such 
disagreement is not found when superradiance physics is 
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P(D) 
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3.4 


1.2 
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5.5 [31, 32] 
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3.3 




18.1 
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4.4 


^^'Na^^Cs 


19.6 
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4.7 


5.2 


33K»5Rb 


26.7 


76 


0.57 [33] 


8.7 


39^133 cs 


30.1 


66 


1.9 f32] 


12 


85Rb"3Cs 


51.8 


49 


1.3 [32] 


21 


Rb (Rydberg) 




33 


~ 10^ 


83 



Table I: Table of the reduced masses n, vibrational spacing 
(jje, dipole moments p, and optical depths OD, for various 
alkali diatomic molecules. To estimate OD , we assume J\f = 
5 X lO'^cm"^ and d = lOO/im. Note that in [H the LiCs 
sample's OD~ 500 is rather small, with a lifetime ~ 20 s. 
The corresponding parameters (n = 40 ^ 39) for rubidium 
Rydberg atoms from [30I ] are also listed for comparison. All 
data are from [sj] (and references therein) except for those 
otherwise specified. 



considered. We are then able to conclude that the same 
cooperativity which leads to superradiance for a higher 
OD still plays an important role here. 

This paper is organized as follows: Sec. [TTldiscusses the 
vibrational transitions in the Dicke-superradiance point 
of view, addressing on the multilevel contribution. Var- 
ious molecules are also compared according to their vi- 
bration relevant parameters in order for better demon- 
stration of superradiance. Sec. IIIII gives the effective 
master equation based on a mean-field approach, where 
the two-body coherence is explicitly kept. Then we show 
our calculation results in Sec. IIVI where we also make 
a comparison for the measured data with our theoretical 
prediction, taking into account the photon and particle 
losses as well as the thermal contribution due to black- 
body radiation. Finally, Sec. |V] summarizes this work 
and future outlook. 



II. VIBRATIONAL TRANSITIONS IN THE 
DICKE LIMIT 

In order to see to what degree the superradiance ef- 
fect must be considered, here we first discuss the optical 
depth (OD) for different setup of molecules. As it is ex- 
pected that a larger OD makes superradiance more ob- 
vious, a closer vibrational spacing is wanted and hence a 
larger wavelength. Therefore those molecules with heav- 
ier reduced mass become better candidates. On the other 
hand, the life time from v to v ~ I transition is usually 
very long. A larger decay rate 7 is preferred so that the 
transition can be observed before other processes such as 
chemical reaction destroy the state. For this reason, a 
stronger dipole moment is preferred. Table U summarizes 
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Figure 1: The enhancement factor {D'^ D~) versus M — J in 
the Dicke picture. Note that the left- (right-) hand side cor- 
responds to upper (lower) levels. Different curves correspond 
to different spin-j cases. Here we have A'^ = 10 and J = Nj. 



the relevant parameters for different molecules which are 
common for cold molecule experiments. Typically the to- 
tal number of molecules is about a few thousands, being 
confined within a trap with a diameter d ~ 10 ~ lOO/iin. 
We also lists the values of OD corresponding to currently 
realizable experiments, showing that quite a range of dif- 
ferent molecules should be expected to have superradi- 
ance. (From this table, one can see that LiCs is not a 
good examplifying system because of its light weight and 
relatively poor OD while KCs and RbCs are best candi- 
dates.) 

Vibrational states form a multilevel ladder where the 
low-lying levels are approximately harmonic. Approach- 
ing the continuum threshold, the level spacing is getting 
narrower when the potential curve becomes more anhar- 
monic. For a typical molecule like LiCs, the number of 
quasi-harmonic vibrational levels can be up to a few tens 
near the potential minimum. In experiments where the 
molecules only have a few vibrational excitations, this 
uniform portion of the vibration ladder should be the 
main stage where spectroscopy is performed. That is, 
most of the vibrational transitions are (nearly) degener- 
ate. In such a situation we are urged to consider multi- 
level effects because the virtual photon exchange, which 
is the main mechanism to build up the Dicke coherence, 
can now take place correlating different levels. 

In order to see the multilevel effect we reexamine the 
enhancement factor {JM\D^ D~\JM) in the Dicke pic- 
ture. Note that in a multilevel case the molecule can 
be described as a spin-j particle. Therefore J — Nj; 
~J < M < J; D+ = Y^i^t with ith particle = 
Y^m=-i 1"^ + The generahzed Dicke states for 



spin-j is given by 



|J,Af) 



' (J + M)! 
{2J)\{J ~M)\ 



{j-y-'^\j,j), (1) 
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where J =J2iJi^ which satisfy 



J±| J, M) = y/J{J +1) - M{M ±1)\J, Mil) 



(2) 



jt \j, m) = + 1) - m{m ± 1) \j,m±l)^ 

with TO indexing the sin gle p article vibrational level. Fur- 
ther calculation yields [27l . [28j the enhancement factor 
{D+D-)jM = N{D+D-) + N{N - 1){D+DJ) with 

{D+ Dr) = l- {j, J -J- -J, M + j|J, ) ' (3) 

and 



{D+DJ)= [{j,J~j;mi,M-rm\J,M) 

mi ,7712 

X (j, J — 2j; TO2, M — TOi — TO2I J — j, M — nil) 

X {j,J - j;mi -1,M -mi + l\J,M) (4) 

X (j, J — 2j; TO2 + 1, M — mi — m2\J — j, M — TOi + 1) 

Figure[T]plots the enhancement factor {D^ D^) jm as a 
function of M — J for various spin-j particles. It is shown 
that the decay rate increases in general as the number 
of levels increases. At M w J, all the curves coincide, 
indicating that they share the same emission behavior 
while starting emitting. The enhancement continues to 
grow as M lowers before it peaks around Af « 0; then it 
turns to decrease with descending M and finally vanishes 
at M = — J. It is worth noting that the curve forms a 
plateau for large j around a range of middle M, and most 
importantly the ascending part of the curve keeps invari- 
ant as j increases. This suggests that the exact number of 
levels is not important in determining the initial emission 
behavior as long as the number of levels is large enough. 
This is valid in typical vibrational spectroscopy within a 
small timescale (compared to 7"^). 

III. TWO-BODY MASTER EQUATION 

The previous approach concerning Dicke superradiance 
in the multilevel case does not take into account the finite 
spatial size of the sample and other decoherence mech- 
anism such as dephasing. These additional factors can 
be contained by considering explicitly the actual dipole- 
dipole interaction. On one hand, dipole-dipole interac- 
tion induced spin (two-level particle) flip-flop, or the so- 
called virtual photon exchange, is the major mechanism 
for generating and maintaining the Dicke type of coher- 
ence. On the other hand, the dipole-dipole interaction 



where r is the inter-particle distance) also intro- 



duces dephasing that suppresses superradiance. These 
effects can be contained in the microscopic Hamiltonian 

2 

H = H^toni + Hj\cM — ^ Pi • Ei(t) — ^pi-Ei(i)), 

i^{l,2} 4=1 



Ho 



where V represents two probe particles [i = 1,2) inter- 
acting with the local field due to environmental particles 
as well as the external field, whose degrees of freedom 
have been absorbed in Hq. In this approach, V is consid- 
ered to be a perturbative term with respect to Hq. Our 
goal is to derive an effective description for the two probe 
atoms in a mean-field manner, which can be thought 
as the next order correction beyond the single particle 
mean-field approximation. Our strategy is to integrate 
out irrelevant degrees of freedom such as field variables 
and environmental dipoles by means of the Schwinger- 
Keldysh Green's function approach. Here a Gaussian 
field is assumed so that we are allowed to only keep up to 
the second order of field correlation. Then the induced 
rate is determined by the field-field correlation. The de- 
tail of derivation can be found in (27i . i28i] and we arrive 
at the effective master equation: 



P = -2^ (pA^ + DiD+p - 2D+pDr 



r + 7 



J2 {pDfD- + D+D-p - 2D-pD+) (6) 



1=1 



- f ^ {pDrD+ + DTD+p - 2D+pD-), 

where p is the two-body density matrix whose dimension 
is (2j -t- 1) X (2j + 1), and 

r = l{e'' - 1)^ + 2C^g'&^Y(t) (7) 



r = 



7^/(C, g) 
r + 7/2 



[?,CQA{t) + 2C'^Q^Y{t)\ 



(8) 



with 



(9) 



(5) 



A{t)= 

m——j-\-l 

i-i 

■^(0 ^ ^ P?n+l,m;m',m' + l ■ 

m,m' = — J 



Here p^^^ = 5 X]i=i 2 tri[p] denotes the reduced single- 
particle density matrix and Pab.cd = c|/9|6, c?) -I- 
{c,a\p\d,h)], where by |a,c) we denote |a)i ® [0)2. For 
both the factor i comes from averaging the interchang- 
ing of two particles. Note that interchange symmetry 
requires pab-cd = Pcd.ab and p*^.^^ = Pba-.dc- The coop- 
erativity parameter is defined as C = 2'KAfc' /loq with 
the particle density AA; g = ujod/{2c) characterizes the 
system size d in units of the radiation wavelength. Func- 
tion /(C, g) = e2«[(C - If + Q^]/{C + Q^?- With these 
definitions, OD = AnCg. 

Each term of master equation ^ can be understood as 
follows: The first line accounts for population pumping 



4 



with the induced rate F. The opposite process, i.e., de- 
cay, is also characterized by F, as presented in the second 
hne. In addition, the free-space spontaneous decay with 
a rate 7 takes place and is contained in the second line. 
Note that at a finite temperature, the thermal contri- 
bution of blackbody radiation may cause extra incoher- 
ent pumping and decay. This effect can be introduced 
by adding an additional rate B « ■f / (^e^'^'^o — 1) with 
= ksT [m to F. Further, the population pumpmg 
corresponds to photon re-absorption (real photons only, 
not virtual photons). This, however, may depends on the 
scattering rate of photons and particles. In order to char- 
acterize the efficiency of re-absorption, we replace in the 
first line F — sF, where the factor s = 1 corresponds 
to maximal re-absorption of emitted photons and s = 
represents loss of emitted photons entirely. Although at 
this point we do not estimate s quantitatively, these two 
extremes s = 0, 1 should be able to give the lower and 
upper bounds for the actual behavior. The third line 
corresponds to generatjon and reduction of inter-particle 
coherence with a rate F. Finally, considering the current 
experimental setup, where there is inevitably particle loss 
out of the trap, one more term —JlP must be added into 
to master equation Note that this term does not pre- 
serve unity trace of the density matrix as time evolves. 



IV. RESULTS AND DISCUSSION 
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A. Multilevel enhancement 

We have investigated the multilevel superradiance in 
great detail in 27|. For completeness, here we include 
a brief discussion for this effect. Such an enhancement 
is due to the fact that the cross-level virtual photon ex- 
change is more possible. In a two-level system, when a 
particle emits a real photon but fails to keep correlation 
to other particles due to certain dephasing, it no longer 
participates in the superradiance coherence unless it gets 
excited later and re-enters the game. This is different in 
the multilevel case because a highly excited particle can 
be involved in manybody correlation for several times of 
its transitions before the ground state is reached. Fig. [2] 
(a) shows the emission intensity /cm as a function of time, 

where /cm = fv-^o J2v '^Pvv corresponding to the overall 
outgoing radiated energy per unit time. These curves are 
dramatically different from an exponential profile yielded 
by single atoms. The initial increase and hump in the in- 
tensity is the signature for superradiance. Note that the 
rising parts of the curves coincide for various level num- 
bers, except that a larger discrepancy is observed in the 
two-level case. This coincidence even extends for longer 
times in the many level cases; the hump is bounded in 
the limit of large level numbers. This can be understood 
by looking at Fig. [2] (b), where we show the top five 
levels' population distribution (diagonal terms of the re- 
duced single particle density matrix) over the vibrational 
ladder at the time when the intensity value is peaked. 



Figure 2: (a) Emission intensity as a function of time. Dif- 
ferent curves correspond to different number of vibrational 
levels. For all data, OD — 3.8 x 10''. (b) Population distribu- 
tion of the five highest levels for 3, 4, 5, and 10-level cases at 
peak intensity (when t = tm). Vm denotes the highest level 
where the whole system is initially populated. Inset: jtm as 
a function of the number of levels. 



We can see that the distribution is almost identical for 
various cases. This happens because the initial emission 
behavior is determined only by the nearby levels. Those 
levels far below do not contribute. 



B. Comparison with experiment 

In usual cold molecule experiments, individual atoms 
are first loaded and cooled in a trap and then through 
photoassociation Feshbach molecules can be formed. In 
most situations, the electronic state of the molecules is 
kept on the ground state but not for vibrational and rota- 
tional states. In a recent experiment done by Deiglmayr 
et al., the vibrational population redistribution of LiCs 
molecules was investigated [3l|. An unusual discrepancy 
was found between the experimental data and the sim- 
ulated curves given by rate equation approaches based 
on Einstein's A-coefficients and inclusion of blackbody 
radiation. Here we apply the superradiance formalism 
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Figure 3: Fitting curves by adjusting the ratio Ni/N-^ be- 
tween the population of levels « = 3 and v — 4. For 
s = 1, i.e., the re-absorption is considered, the optimized ratio 
N4/N3 — 7.47; for s = 0, i.e. the re-absorption is neglected, 
the ratio N4/N3 = 6.72. An optimization is also taken by 
overall rescaling the y value for best fit. For comparison, 
two simulated curves from [sj are plotted, corresponding to 
inclusion of the first-order and second-order particle losses, 
respectively. The thin-solid line corresponds to the case when 
the balckbody radiation is considered, with the optimized ra- 
tio N4/N3 — 7.55. The resultant difference, however, is rather 
small. The scattered dots with error bars are the experimen- 
tal data. 



mining the relaxation behavior. It is notable that, in this 
experiment, OD « 500 is very low, even lower than the 
values suggested in TableU] It is actually near the bound- 
ary of superradiance and ASE. Note that ASE as well as 
superradiance are both caused by the dipole-dipole in- 
teraction. But in ASE the decoherence dominates and 
therefore the inter-particle coherence is less pronounced. 
It should be emphasized that the ASE mechanism is also 
contained under the framework of our approach. In fact, 
LiCs is a rather bad choice for demonstrating superradi- 
ance because of its low reduced mass. A heavier species 
such as RbCs under similar experimental setup will do 
better (see Table |l|. 



C. Blackbody radiation 

An additional curve (thin-solid line) is also plotted in 
Fig. [21 taking into account the thermal contribution of 
blackbody radiation. But the deviation is rather small. 
Here we briefly discuss the role of the blackbody radiation 
repumping and decay. In Eq. ^ it serves as a factor in 
addition to F. At room temperature T — 300K, we can 
estimate B = O.757. This is usually a significant quantity 
in the single particle case but not in the superradiance 
case, where F is typically 2 orders of magnitude larger 
than 7 for OD = 10^ ~ lO^. 



to their experimental setup, where OD « 500. In Fig. 
13] we show the time evolution of the population profile 
corresponding to level v = 3, directly comparable to the 
photoassociation measurement outcome. 

We fit the experimental data as follows: From [3l| we 
find the ratios of the lowest ten vibrational levels. In 
the experiment, however, most electrons (^^70%) are still 
in some of the upper levels, with unknown distribution. 
Thus, we have to allow for at least one fitting parameter 
to take this unknown into account. In our case, for ease 
of calculation, this fitting parameter is the ratio between 
the population of levels v — 3 and w = 4. The curves are 
then calculated for 51 total levels and an approximation 
where only the lowest 10 of those 51 levels are used (Fig. 
[3] shows that there is nearly no difference between those 
two cases.) 

In Fig. [31 the scattered dots with error bars repre- 
sent the experimental data and the thick-solid (thick- 
dashed) lines correspond to the overall 51-level (10-level 
truncated) case. For comparison, we plot s — and 
s = 1 cases corresponding to the upper and lower bounds, 
respectively, when re-absorption efhciency is or maxi- 
mized. We find that the s — curves show a very good 
agreement with the experimental data, suggesting that 
re-absorption is not a significant effect when the trap is 
transparent to radiation. Most importantly, such agree- 
ment implies that superradiance (or at least the inter- 
particle cooperativity) plays an important role in deter- 



V. CONCLUSION 

The vibrational relaxation behavior of cold molecules 
displays superradiant behavior when their density is large 
so their particle cooperativity must be considered. We 
have developed a mean-field method that explicitly take 
into account the actual dipole contribution as well as 
the multilevel effect. Our theory works rather well de- 
scribing the experimental data from measuring the relax- 
ation of vibrational population for cold molecules. This 
indicates that the cooperativity has a significant effect 
regarding spectroscopy in most ultracold molecules ex- 
periments, including vibrational and/or rotational tran- 
sitions. Our calculations suggest that considerable co- 
operative effects can be expected from a large amount 
of ultracold molecule experiments if vibrational and ro- 
tational transitions are considered. In general, for all 
setups with an OD density of the order of 10^ ~ 10'^ or 
higher, superradiance will play an important role. 

Some open questions remain such as how superradi- 
ance is modified by (1) the anharmonicity of the vibra- 
tional ladder, (2) the non-constant transition dipole mo- 
ments from level to level, and (3) the trap-induced inho- 
mogeneity in density. They are, however, at this point 
not expected to have a major qualitative impact on the 
results. 
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